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Abstract 

We present a generalization of the Hartree-Fock Bogoliubov (HFB) theory in which the coupling 
between one and two quasi-particles is taken into account. This is done by writing the excitation 
operators as linear combinations of one and two HFB quasi-particles. The excitation energies 
and the quasi-particle amplitudes are given by generalized Bogoliubov equations. The excitation 
spectrum has two branches. The first one is a discrete branch which is gapless and has a phonon 
character at large wave-length and, contrarily to HFB, is always stable. This branch is detached 
from a continuum branch whose threshold, at fixed total momentum, coincides with the two quasi- 
particle threshold of the HFB theory. The gap between these two branches at P=0 is equal to 
two times the HFB gap, which then provides for the relevant energy scale. We also give numerical 
results for a specific case. 
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INTRODUCTION 



The experimental realization of Bose-Einstein condensation in trapped neutral bosonic 
atoms has opened the opportunity for a comparison of microscopic theories of dilute systems 
with experimental data The standard approach is to solve the Gross-Pitaevskii (GP) 
equation for the condensate wave function and the linear Bosroliubov-de Gennes (BdG) 

Q 

equations for the collective excitations For T = this theory has been successfully 
compared with the existing data 

Physically, the GP + BdG theory is a free quasi-particles theory. However it is not the best 
of such theory from a variational point of view, being superseded by the Hartree-Fock- Bo- 

n 

goliubov (HFB) theory where the non-linear BdG equations are solved self-consistently 
and the GP and BdG equations are coupled. The two approaches agree when we neglect 
the depletion of the condensate in the self-consistent theory. This would suggest that 
HFB is the proper theory to use when the fluctuations become important, for example, in 
trapped gases through the mechanism of Feshbach resonance or at finite temperature. 
It is however well known that the HFB theory has problems when applied to homogeneous 



le Hugenholtz-Pines 
. The HFB theory 



systems p . Indeed, the excitation spectrum has a gap which violates t 
theorem, which states that the excitation spectrum should be gapless 
is also in contradiction with some thermodynamic data in ^He (for example, the specific 
heat), which requires that, in the large wave-length limit, the excitation spectrum should 
have a phonon behavior j^. 

In this paper we present a generalization of the HFB theory leading to an excitation 
spectrum which, in particular, eliminates the gap problem. The key ingredient for this is the 
inclusion of two quasi-particle components together with their coupling to one quasi-particle 
components for the description of the excitation spectrum. Excitation energies and the 
structure of the corresponding modes are given by generalized Bogoliubov equations. The 
excitation spectrum includes a discrete stable phonon-like excitation branch. In addition 
to this discrete branch there is a continuum branch whose threshold is located at twice the 
HFB quasi-particle energy. This branch has, therefore, a gap at total momentum P = 
which is twice the HFB gap. Therefore, this gap sets an important energy scale for the 
excitation spectrum. 



2 



The use of this mechanism to solve the HFB gap problem has in fact being pioneerec 
more than for decades ago by Takano jl^. A more recent work by Hutchinson at al 
deals with the coupling between one and two quasi-particle components in a perturbative 
way, while the work of Kerman and Tommasini ^U] deals with the same problem on the 
basis of the Gaussian functional approximation to a field theoretical variational procedure. 



In this paper we use standard equations of motion techniques [13|, [IJ], employed, for 
example, in reference to study chiral symmetry restoration in the linear a model 
non-pertubatively. This allows for a clear identification of the dynamical role played by 
various parts of the many-body hamiltonian when expanded in terms of quasi-particles. 



Our paper is organized as follows: In section 2 we briefly discuss the basic properties of the 
HFB theory. In section 3 we derive the generalized Bogoliubov equations for the excitation 
energies and quasi-particle amplitudes. Our derivation allows for a clear identification of 
parts of the hamiltonian responsible for the coupling between on and two quasi-particle 
components of the excitation operator. In section 4, we show that there exist a Goldstone 
mode at momentum equal to zero. Our proof of its existence is very simple and clearly related 
to the violation of number conservation. The results of a numerical application of the theory 
are discussed in section 5. Specifically we examine the properties of the excitation spectrum, 
its stability, and change in physical content as a function of the total momentum P. We 
also make a comparison with the HFB and Bogoliubov approximations. Our conclusions 
are presented in section 6. All the expressions needed for numerical applications are given 
in appendices A and B. 

THE HFB THEORY 

The starting point is the Grand-Hamiltonian written in second quantization as 

h = H - fiN = J2{ek - fi)ala^ + ^ ^ ^(?)4i+q42-q«kiak2 (1) 

k ki,k2,q 
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where is the free particle kinetic energy, 

V{q) is the Fourier transform per unit volume of the atom-atom interaction potential 

V{q) = i / V{r)e^^-^£r = (3) 

and the operators aj^ e a^, respectively, create and annihilate atoms in a state with mo- 
mentum hk, the corresponding wave-function exp(2k.r)/v^ satisfying periodic boundary 
conditions in a cubic box of volume F. 

In a first step we perform a canonical transformation to the quasi-particles by introducing 
a new set of creation and annihilation operators through the Bogoliubov rotation 

flk = Ck + ^o5k,o = ^ik?7k - VkV-k + zoSkfi, (4) 

where and fk are even functions of k, Uk = u_k, fk = "^-k; and Zq is a c- number. The 
constant Zq appears as a shift in the equation for k = to account for the macroscopic 
condensate in the zero momentum state. In order to render the transformation canonical, 
the Bogoliubov factors have to obey the constraint 

ul-vi = 1. (5) 

It is straightforward to write the Grand-Hamiltonian h in the quasi-particle basis. After 
normal ordering one obtains, 

h = H - nN = hQ + hi + h2 + h^ + K (6) 

where the normal ordered operators hi contain i quasi-particles. They are written explicitly 
in Appendix A. 



The amplitudes tik, Wk and the shift zq are determined in the HFB theory by minimizing 
the expectation value of /i, ($|/;.|<l') in the quasi-particle vacuum, that is, ?7k|$) = 0. Since 
the only term which contributes to the expectation value is the term /iq, the minimization 
is equivalent to the equations 
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dzo 



(7) 



dho dhoduu 



with ho given by 



ho = -ZoV + 1^(0) {V{0) + V{k))z',]vl 

1 



2 

k ki,k2 
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+^ E ^(1^1 - k2|)Mkit;ki^ik2^^k2 (9) 

ki,k2 

The two equations above can be written in a very compact way if we introduce the 
Hartree, exchange and pair potentials defined by the following relations jl?! : 



= E^(o)(<fK«k.l$) (10) 

ki 

UM = En|k-ki|)($|at^ak,|$) (11) 

ki 

Upairik) = E^(|k-ki|)($|ak,a_k,|$) (12) 

ki 

These potentials can be written as the sum of two terms which can be interpreted as related 
to the condensate and to the non-condensate, respectively 

U, = [/,^ + t/r = nO)^o' + EnO)($l4ck,|$) (13) 

ki 

[/e.(k) = [/:, + ^- = V(A:)^o' + En|k-ki|)($|4ck,|$) (14) 

ki 

U,a^r{k) = f/;,.. + f/,":,. = ^^o' + E ^ (I k " kl |) ($ | Ck, C_k J $) (15) 

ki 

In terms of these potentials the equilibrium equations (|7j) and (jH)) can be written as 

zo[-fi + u, + u:m + u;:,M] = (16) 

tanh2a, = ^-%i\^ = (17) 

ek + Uh + Uex[k) - /i e(k) 
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with Mk = coshcTk and fk = sinhcik. 



The quasi-particle vacuum |$) does not have a definite number of particles. In order to 
control the number of particles we determine fi from the condition that the mean value of 
the number of particles in the state |$) is A^, = ($|iV|$) which gives 

N = zl + Y.vl (18) 

k 

Thus, the set of equations (fTBjl . (fTTj) and (fTH|) determine /i, Zq and the Bogoliubov 
amplitudes and v^. Equations ( ITBj) and (fTTj) can also be derived by demanding that hi 
vanishes and that /i2 is diagonal in the quasi-particle basis, /12 = J2k^(M)VkVii^ with ti;(k) 
the quasi-particle energies, ci;(k)^ = e(k)^ — UpairiX)^- 

One feature of the HFB theory is that the excitation energy have a gap in the limit k —>■ 

B. 

coiof = -W{0)noU;:„{0) (19) 

where hq = Zq /T is the condensate density. The existence of an energy gap in the excitation 
spectrum does not agree with the phonon spectrum in superfiuid systems and is also 
in contradiction with the Hugenholtz-Pines (HP) theorem, which sates that the en- 
ergy of an excitation with wavenumber k of a many-body system should vanish when /c — 0. 



An approximate way to satisfy the HP theorem is to neglect the so-called anomalous 
density rrik = ($|ckC_k|'^') in the HFB theory [P]. In this approximation Up^^^(k) vanishes 
and the gap disappears. This approximation is known as the Popov approximation. In the 
next section we are going to present a theory that has a gapless dispersion relation while 
taking mk fully into account. As it turns out, this theory also gives a physical meaning to 

(uni). 



THE QUASI-PARTICLE RPA 



As is well known in many-body physics the Random Phase Approximation (RPA 



out the Goldstone mode due to a symmetry breaking at the mean field level (l3, 



li 



singles 



One among the many ways of deriving the RPA equations, is the hnearization of the 
equations of motion [2^. In principle, if we find operators satisfying the equations 



[H, Qi] = n^Qi (20) 

Qal^o) = (21) 



with the normahzation condition 



{^o\[Q^,Qimo) = Sfs,a (22) 

where |\l/o) is the exact ground state, we have found an exact excited state of the many-body 
system since, from the above equations, it follows that |\1/q,) = Ql^l'^o) is an eigenstate of 
H with excitation energy equal to Qa- However this cannot be carried out in general, and 
we are bound to use approximate methods regarding both Qa and |\I^o) in the solution of 
equations ()20j) and (PT|) . In the method of the linearization of the equations of motion we 
make an ansatz about the form of the excitation operators, writing it as a linear combination 
of basic excitations and we linearize the left -hand side of equation ()20|) with respect to these 
operators. 

In this paper we look for excitation operators which are a combination of one and two 
HFB quasi-particles, 

Vp — XpVp + VpV-P + ^q,P 1 r -'q.P / (.^-jj 

q>0 vl+^q,0 Vl + <^q,0 



In the equation Qp creates an excitation with momentum P and it is a linear 
combination of one, r^p, ?7_p , and two, ''7q+p/2^-q+p/2' ''?-q-P/2'7q-p/2 HFB quasi- particles. 
This last term creates (annihilates) a pair with total momentum P(— P) and relative 
momentum q. The coefficients xp, y-p, ^q,p and Yqp are even functions of P. As 
the pair creation and annihilation are invariant by the replacement q —>■ — q, Xqp and 
Yq^P are even functions of q and we restrict the sum in order that the pairs appear only once. 

The coefficients in fl23D are determined by the method of the equations of motion in the 
version of references jlS^and jl^ . which is a systematic way of achieving the linearization 
referred above. Following these references, notice that from equations (|^m) and (pT|) one has 



(v|/o|[Qp,i/,Ql,]|vI/o) = Qp{^o\[Qp,Qlm (24) 

where [A, is the symmetrized double commutator ^{[A, [B,C]] + [[A, B],C]). 
Requiring that Qp be stationary in a variation of the excitation operators one has: 

{nsQp,H,Qim = npm6Q^,QU\<^) (25) 

where 6Qp is the hermitian conjugate of the operator given by the variation of the coefficients 
in and, as usual, we replaced the ground state |\E'o) by the HFB vacuum , |$). 

Performing the variation indicated in equation (j2Hl), we get for the excitation energies Qp 
and the coefficients 

X =\ \ and y= , (26) 



p 




the equation 



B* A* ) \ y J \o -1 ) \ y 

with the coefficients obeying the normalization condition 



(27) 



($1 [Q^, Q^^) = x^'x^ - y^'y^^ + K'p^^p " ^S^n,p) = '^A.r- (28) 

q>0 

For each P we have a number of modes equal to the number of operator pairs plus one, 
n-pairs + 1, actually this number is denumerable infinite, and this is denoted by the quantum 
numbers A, r in equation (j2HI)- 

In expression (j^ . Xp(Fp) stands for the set of Upairs coefficients Xqp(Fqp) and A and 
B are, respectively, hermitian and symmetric matrices of dimension Upairs + 1- 

The hermitian matrix A is given by 

^= I '1 '1 I- (29) 




The diagonal blocks A^^ and A'^'^ are hermitian matrices with dimensionality 1 and np^ 
respectively. They are given by 
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A'\p)^{nv^,H,nim, (30) 

and 

A''{c^, q; P) = H, !^Zi2k±^|$). (31) 

The coupling matrices A^^ and A^^ are hermitian conjugates with dimensions 1 x ripairs and 
^poir-s X 1 respectively and are given by 

^-(q;P) = ($|[,p,if,%a2k^]|$), (32) 

V 1 + <^q.o 

The symmetric matrix B can be split in a similar way 

(33) 



The diagonal blocks B^^ and B'^^ are symmetric matrices whose elements are given by 



yB^' B^^ J 



B'\P)^mrii>,H,V-^m, (34) 

and 



The coupling matrices B^^ and are transpose of each other and given by 

B^^(q;P) = ($|[77p,if, ^-^-^/^y^/^ ]|$), (36) 

+ V 

Note that the matrices 12 and 21 couple one and two quasi-particle excitations whereas 
the matrices 11 and 22 act only inside the one and two quasi-particle basis, respectively. As 
shown in appendix B, the matrices 12 and 21 depend only on h^, (see Eq. 6). Therefore this 
term is responsible for the coupling between the one and two quasi-particle components. 
On the other hand the matrices 11 and 22, that do not mix these two components, depend 
only on ^2 and h^. All the matrix elements are computed in detail in appendix B. 
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If the coupling terms /13 and /i4 are set to zero we have that B = and A is diagonal 
with eigenvalues uj(P) and co'2(q, P) = i-i;(q + P/2) + ti;(— q + P/2) which correspond to one 
and two free quasi-particle energies. 

As the static quantities are real, the elements of the matrices A and B are real and the 
RPA equations can be written in a more compact and symmetrical form if we introduce the 
following new variables: 



MP) 

02(q;P) 
vri(P) 
vr2(q;P) 



(ATq^p — V'q^p) 



(37) 



Grouping 0i,02 and tti, 712 in order to have the following ^^p^jj^g + 1 column matrices 



Q 



' 0i(P) \ 






v</'2(q;P) ] 





^l(p) 

^2(q;P)J 

we can rewrite the equations in a very compact form jlll . Vi 



(38) 



= AJP (39) 

= A^Q (40) 

where now we work with two symmetric matrices ,4+ and A- which are given in terms of 
the matrices A and B as 

A+ = A^B (41) 

A- = A-B (42) 

The elements of these matrices can be written as: 



A^l = A 



11 



(43) 
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^i'(q, P) = , '\ go(P){Co(-)Co(+)[l^(+) + V^(-)] + [Co(+)Co(-) 



1 + 5q,0 

-C^\+)C^\-)]Vm (44) 

^+(q,P) = , '\ c,\v){[c,\+)c,{-) + Co{+)c,\-)][v{+) + 

2W1 + dq,o 



^(-)]} (45) 
^?.2(q,q',P) = [^(+) + ^(-)]5q,q,+ 



2^(l + M(l + V,o) 
+C^o(+')Co(-')Co(+)Co(-)][F(|q - q'l) + \/(|q + q')] + [Co-^(+')Co-^(-') - 

C,{V)C,{-')][C,\+)C,\-) - Co(+)Co(-)]l-(P)} (46) 

^^2(q,q',P) = [cu(+) + cu(-)]5q,q,+ 

:[Co-^(+')Co(-')Co-^(+)Co(-) + 



2^(l + 5q,o)(l + V,o) 
Co(+')Co-^(-')Co(+)Co-^(-)]r(|q - q'l) + \C,\+')C,{-')C,\-)C,{+) + 

Co(+')Co-^(-')Co(-)Co-^(+)]V(|q + q'|). (47) 

These expressions show the remarkable result that all the matrix elements depend on the 
static factors through only one function i.e. 

Co(q)=Wq-t^q (48) 

where we used the notation ± = ±q + P/2 and ±' = ±q' + P/2. By solving the coupled 
equations (j39|) and (j4m) we will find the new excitations which will now have one and two 
quasi-particle contributions. The results for the excitation energies will lead to a discrete 
branch and a continuum whose threshold coincides with the two quasi-particle threshold 
of the HFB theory. The discrete branch is detached from the continuum and due to the 
coupling between the one and two quasi-particles will be pushed down and become gapless. 
This fact can be proved for any pseudo-potential as will be shown in the next section 

THE GOLDSTONE MODE 

Equations (fTTj) have a zero energy solution when P = 0, fio = 0, with zero-norm, 
the Goldstone mode To identify this solution one has to consider the generator of 
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the symmetry violated by the theory which, in our case, is the U{1) symmetry whose 
generator is the number operator. This operator when written in the HFB basis possesses 
components that are present in the general "anzatz" for the excitation operators, eq.(j22I)- 
These components will be identified with the excitation operator of the Goldstone mode Qg- 

Thus to find Qg we start writing the number operator 

^ = E4«k (49) 

k 

in terms of the HFB quasi-particles r]l., r^k, giving 



N = N + zoiuo - vo){vo + ^o) - E rTT^(^I^-k + ^k^^-k) + E(«k + ^Dvlvi.- (50) 

k>o J- + ^k,0 k 

Comparing with the general "anzatz" eq. ()23|) we identify the excitation operator of the 
Goldstone mode with 

Qg = Zo{uo - Vo){r]o + Vo) - (^k^-k + ^k^-k)- (51) 

k>o ^ + "^k.O 

Given this form, our next task is to prove that, when P = 0, there is a zero-energy solution 
of ()27|) of zero norm with. 



xo = yo = zo{uo - vo) (52) 

-'^q.O = ^q,0 = 7===- (^^) 

y 1 + Oq.O 

Equation is equivalent to the coupled equations (jH^ and (jin|) . Since the Goldstone 
mode has zero-energy and zero-norm one has fio = and V = and the coupled equations 
reduce to 



A+Qg 



(54) 



with 



Q 



G 




2xo 



(55) 
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From the expression of the matrix elements of ^4. , Eqs. (44) and (46) (at P = 0), Qg given 
above is a zero-energy solution of the eqs. (fH^ and pUj) provided the following identities 
are satisfied 



lu,{uo - voY + U;:,M = (56) 

Upair {(l) - 2t^qUqt;q = (57) 

These identities are easily seen to hold if we use the following relations obeyed by the static 
quantities: 



<-\[^,^')' '''' 

2u^v^ = ^Ip!^. (59) 

In the equations of motion method the connection between the Goldstone mode excitation 
operator Qg and the number operator N goes as follows. Since [h, N] = one has 

msQ,[h,N]m = o. (60) 

At first glance there is a difficulty to conclude from the above equation that Qg is a zero- 
energy solution of Eq.(j2Sl), caused by the presence of the term I]k(^k + '"kJvl.Vk Eq. ljMH) 
which does not belong to the general "ansatz" eq. ()23|l . However this term does not give 
any contribution to ()6()j) and since in our case the double-commutator is identical to the 
symmetrized double-commutator one has 

msQ,h,QGm = o (61) 

showing that the last term in Eq.(48) does not play any role and indeed Qg is a zero-energy 
solution of eq.(21). 

In the HFB case we could proceed in the same fashion. However in this case the terms 
which do not belong to the HFB "ansatz" do contribute to the matrix element (25) and, as 
a consequence, the HFB equations do not have a zero-energy mode, the excitation spectrum 
always has a " gap" . 
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NUMERICAL RESULTS 



In this section, to ilustrate the predictions of the theory outhned in this paper and 
compare with the HFB and Bogohubov approximations, we present and discuss the results 
obtained by solving equations pUj) and HOI in a specific case. 



as our ps 

HQ 



To begin with, we choose as o ur p seudo-potential in momentum space, V{q), a purely 
repulsive Gaussian defined by 



V(k) = e-— (62) 

m 

where, as usual, the pseudo potential at = and the scattering length are related by the 
expression 

mVjO) 

" = ^''^ 
We measure energies in units of /i^/(2ma^) and lengths in units of the scattering length 
a. The width of the pseudo potential was chosen to be of the order of the scattering length 
a = 2.8a as in [2^. The only parameter left is the total density p which was taken to be 
such that a^p = 10~^ and 10~^. These choices of the "dilution parameter" falls between 
the values corresponding to the dilute experiments and to liquid Helium. Values such as 
10^^ and 10^^ for the trapped EEC can be achieved in experiments conducted close to a 
Feshbach resonance 

Once the parameters are specified, we calculate the energy of the discrete branch, the 
continuum threshold and the structure of the excitation operators of the discrete branch. 
The first step in these calculations is to solve the self-consistent static equations, (fT^ . lfTTj) 
and ()18|) which are needed in the determination of the and A- matrices ()41|) and ()42j).The 
next step is to solve the coupled equations (jniJ|) and (jlOj) . The standard way to proceed is 
to take the thermodynamic limit and solve the corresponding coupled integral equations. 
In this paper, we took a different route: we solved the coupled equations (jHEjland pHjl .for 
a box with volume V. The value of the volume V is increased and the whole calculation 
is done again, until "saturation" is observed, indicating that the thermodynamic limit has 
been sufficiently reached for these quantities. 
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In Fig.l we present results for a?p = 10~^ since the qualitative behavior of the calculated 
quantities does not depend on the value of the density. We start by looking at the discrete 
branch in Fig. la. In the long- wave length limit this branch is gapless, as shown in section 
3, has a phonon like dispersion relation i.e Qp = cP and it is always stable. 

We have also verified numerically that the continuum threshold starts at the minimum 
value for the two quasi-particles HFB energies, t^2(q, P) = c<j(q + P/2) + uj{—q + P/2) at a 
fixed value of P, which in our case always happens at q = 0. 

We can compare our results with the HFB and Bogoliubov approximations. In the HFB 
approximation we have free quasi-particles, and the one and two quasi-particle branches are 
decoupled. As shown if Fig. 1(b) both branches have a gap and, in the limit of long wave 
length, it is not linear in P. It is possible to show that these two branches always cross at 
some value of P. At small P the one quasi-particle branch energy is always lower than the 
two quasi-particle threshold due to the existence of the gap, 002(0, 0) = 2uj{0), whereas for 
large P is just the opposite, co'2(0, P) ~ < uj(P) ~ P^, always leading to a crossing of 
the one quasi-particle branch and the two quasi-particle threshold. As a consequence of the 
crossing the one quasi-particle branch always becomes unstable eventually. 



For a^p = 10"^, the crossing point happens at P = .63a^^ therefore the one quasi-particle 
branch is stable for P < .63a~^, becoming unstable for P > .63a~^. If we compare with Fig. 
la we see that the discrete branch "avoids" the crossing moving away from that point and 
rapidly approaching the two quasi-particle threshold after the crossing point. This effect is 
seen in grater detail in Fig. 2(c). 

In the Bogoliubov approximation, shown in Fig. 1(c), the two branches are gapless 
and phonon like in the lon g w avelength limit. In this approximation the one quasi-particle 
branch is always unstable 2l|.For a^p = 10~^ and momenta P < 0.54a~^ the one quasi- 
particle branch and the two quasi-particle threshold are degenerate. In this case the one 
quasi-particle decays into two quasi-particles where one of them carries all the momentum 
and energy. This is possible because the one quasi-particle branch is gapless. For P > .54a~^ 
the one quasi-particle branch lies above the continuum threshold that occurs for zero relative 
momentum q = 0. 

In conclusion, we found that the energy of the discrete branch, interpolates between the 
Bogoliubov one quasi-particle spectrum and the HFB two quasi-particle threshold, with a 
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FIG. 1: The excitation spectrum according 
to the RPA, (a), to the HFB theory, (b), 
and in the BogoUubov approximation, (c) 
for a^p = 10~^. In fig 1(a) the sohd lines 
indicate the discrete branch and the contin- 
uum threshold, whereas in figs (b) and (c) 
the lower solid line is the one quasi-particle 
branch. The dashed line corresponds to the 
region where the energy of the one quasi- 
particle branch is greater than the thresh- 
old of the two quasi-particle branch. The 
excitation energy is measured in units of 
h'^ /{2ma?) and the momenta in units of h/a. 
See text for details. 



a p=10 




relatively sharp transition region near the onset of instability of the HFB one quasi-particle 
spectrum, as illustrated in Figs. 2(a), 2(c) for o? p = 10~^ and in Figs. 3(a), 3(c) for 
a^p = 10^^. From these graphs we also see that the sound velocities are practically equal to 
the sound velocities of the Bogoliuvov approximation. 

We extracted information on the composition of the excitation operator as a function of 
the total momentum P by calculating the quantity 

c,CP)=xl-yl, (64) 

which corresponds to the relative weight of the one quasi-particle component. Analogously 
for the two quasi-particle component we can define 



C2(P) = E«P->^,P)- 

q>0 



(65) 



These two relative weights are related through the normalization condition ()28|1 that gives 



Ci(P)+C2(P) = l 



(66) 



Fig. 2(b) and 3(b) show the values of Ci(P) and C2(P) as a function of the total momentum 
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FIG. 2: Fig. 2 (a) shows the RPA discrete 
excitation energy and continuum threshold 
(soUd hues) and the Bogohubov one quasi- 
particle excitation energy (dot-dashed hne). 
In Fig. 2(b) ci(P) (sohd hne) and C2(P) 
(dashed hne) give , respectively, a measure 
of the one and two quasi-particle character 
of the discrete branch, as a function of P, for 
a^p = 10~^, Fig. 2(c) is a zoom of Fig. 2(a) 
near the transition region. Units as in Fig.l. 



P for a^p = 10^^ and 10^3 respectively. In the long wavelength regime the excitation 
operator is predominantly a one quasi-particle operator. On the other hand in the short 
wavelength regime it becomes predominantly a two quasi-particle operator as it approaches 
asymptotically the continuum threshold. Note that there is a sharp transition between 
this two regimes, the effect being more pronounced at higher densities. Comparing Figs. 
2(a), 2(b) and 3(a), 3(b), we see that the change from a predominantly one quasi-particle to 
a predominantly two quasi-particle physical content of the discrete branch of the excitation 
spectrum and of its corresponding excitation operator occur in the same momentum range. 

CONCLUSION 

In this paper we presented a generalization of the HFB theory, cast in terms of well 
known methods of equations of motion in order to access the excitation spectrum of a 
condensed many- boson system. The key ingredient for the generalization is the coupling 
of the one and two quasi-particle components to form the excitation modes. These are 
determined by solving the appropriate generalized Bogoliubov equations for the relevant 
amplitudes and excitation energies. 
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FIG. 3: Fig. 3 (a) shows the RPA discrete 
excitation energy and continuum threshold 
(soUd hues) and the Bogohubov one quasi- 
particle excitation energy (dot-dashed hne). 
In Fig. 3(b) ci(P) (sohd hne) and C2(P) 
(dashed hne) give , respectively, a measure 
of the one and two quasi-particle character 
of the discrete branch, as a function of P, for 
a^p = 10'^, Fig. 3(c) is a zoom of Fig. 2(a) 
near the transition region. Units as in Fig.l. 



The generalized Bogoliubov equations are shown to have a Goldstone mode at zero 
transfered momentum, whose structure is related to that generated by the particle number 
operator. Correspondingly, an examination of the properties of the excitation spectrum 
reveals a detached gapless excitation branch with phonon-like dispersion at small momenta 
and a continuum branch starting at the two quasi-particle threshold of the HFB theory. 
The detached branch has a predominantly one quasi-particle character at small momenta, 
where it closely approximates the Bogoliubov one quasi-particle spectrum. At high 
momenta it approaches the continuum threshold and eventually acquires a predominantly 
two quasi-particle character in what can be seen as an avoided crossing situation, due to 
the coupling between one and two quasi-particle components included in the calculation. 
The transition from one to two quasi-particle character is relatively sharp and occurs 
near the onset of instability of the HFB one quasi-particle spectrum. Differently from the 
Bogoliubov and HFB approximations the detached phonon-like branch is always stable. 

The presence of the continuum branch serves, in particular, to give physical significance 
to the HFB one quasi-particle energy gap, this being the quantity which sets the appropriate 
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energy scale for the continuum threshold. 

The features revealed in this generahzed theory are closed hnked to the depletion of 
the condensate caused by the two body interaction effects. They should be particularly 
relevant, therefore, in cases where such depletion effect become important. This happens for 
larger condensates densities and/or larger values of the relevant scattering length (implying 
stronger effective interaction) a situation that may be realized experimentally, for example, 
by taking advantages of Feshbach resonances. 
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Appendix A 



THE GRAND-HAMILTONIAN WRITTEN IN NORMAL ORDER 

The Grand-Hamiltonian can be written in the following form 

h = ho + hi + h2 + h + hi (67) 

where hi corresponds to the normal ordered component with i quasi-particle operators. In 
turn,we will write each hi as : 

i 

j=0 

where i — j stands for the number of quasi-particle creation operators and j for the num- 
ber of annihilation operators. In what follows we will give explicit expressions for all the 
components appearing in (jU7|) 

i) ho 



ii) hi 



where 



ho = -zl^i + $V{0)+Y.[ek-^i+iViO) + V{k))z',]vl 
^ k 

k ki,k2 

+ ^ E ^(|kl - k2|)Mki^^kiMk2^k2 (69) 
ki,k2 



h\ = hiQiTjl + hQir]o (70) 



/^lo = V = (-/i + Uk + t/:;(o) + u;:,MU{xo - y^) (7i) 

From eq. (f7T|) it follows that /ii = as a consequence of the equilibrium equation (fT^ . 
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iii) h2 



h = ^/iii(k)r7[r7k 

k 

+/i2o(k)r7lr/lk + /io2(k)?7kr/-k (72) 

where 

/iii(k) = e(k) cosh2crk — Upair(k.) sinh2(Tk (73) 

/i2o(k) = /io2(k) = -e(k) sinh2(Tk + ?7pair(k) cosh2o-k (74) 

From the equihbrium equation (fT7j) we see immediately that the non-diagonal components 
in equation (ff^ . /i2o(k), /io2(k) vanish and the coefficient of the diagonal component /iii(k) 
is equal to cuk leading to 

/i2 = E^(k)^I^k (75) 

k 

iv) /l3 

/is can be written as 

^3 = X!{^i.2(k,k')r?k+k'7k^?k' + /i2,i(k,k')?7k/?7kr/k+k' 

k,k' 

+^3,0 (k, k')r?I+k'^-k^lk' + ^0,3 (k, k')r/k+k'?7-k?7-k'} (76) 
where the hi k coefficients are given by 



/ii,2(k, k') = /i2,i(k, k') = y {K+k'Mk' + ^^k+k'^^k'] Nk - vi,]V{k) 
+ K+k'Mk + ^^k+k'^^k] Nk' - vu']V{k') 

-K+k' - ^^k+k']Nkfk' + fkMk']V^(|k + k'l)} (77) 
/i3,o(k, k') = /io,3(k, k') = y {[wk+k'fkfk' - ^^k+k'MkMk'][V^(/i;) + V{k')]} . 

Note that all the coefficients obey the symmetry propriety hij(k,k.') = k, — k') = 
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v) h4 

Analogously to the previous case we can write as 



h4 ^ ^I,3(kl,k2,q)77j^2_q77-ki-q77ki?7k2 
ki,k2,q 

+/i3,i(ki, ka, q)r^ki^L^-ki-q^k2-q 
+/i4,o(ki, ka, q)?7ki+q?7k2-q^-ki?7lk2 
+/io,4(ki, ka, q)r/ki+q??k2-q^-ki??-k2 
+/i2,2(ki, ka, q)r/ki+qr/I,_qr/kir;k2 

where the coefficients are given by 

^3,i(ki,ka,q) = /ii,3(ki, ka, q) = -Ka-q^'ki+qMkiMka + 

^^ka-q^iki+q^^ki-yka] ^(g) 

^4,o(ki, ka, q) = ^o,4(ki, ka, q) = ^Uk^+qUk^-ciVk^Vk^V (q) 

1 

/i2,2(ki,ka,q) = 2{[^ki+qMki + ^'ki+q^'kill^ka-q^ka + ^^ka-q^^ka 
+ [lik2-qtiki^^ki+q^^k2 + ^ikl+q^ik2^^k2-q^^kl]^(|kl + ka|)} 

These coefficients,except h^^i and /ii,3 obey the symmetry property 

hij{ki, ka, q) = /iij(ka, ki, -q) = hij{-ki, -ka, -q) 
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Appendix B 



RPA MATRICES 

In this appendix we evaluate the elements of the matrices A and B. The key property 
that we will use to calculate the average value of the symmetrized double commutators is 
that 1$) is the quasi-particle vacuum. In order to work with more compact expressions we 
introduce the following quantities 



where (|H3j) and (jHUj) are usually called coherence factors and are well known in the study of 
Bose systems 

The matrix A 

^^^(P) is given by, 



Ci(k, k') = Ci(k', k) = Cii-k, -k') = MkWk' + v^v^, 

C2(k, k') = C2(k', k) = C2(-k, -k') = Mkt^k' + t^kUk'. 



Co(k) = Mk-Wk 



(84) 
(85) 
(86) 



A'\P) = mVP,H,rjl,]\<^). 



(87) 



From the property that |$) is the quasi-particle vacuum it follows immediately that only /i2 
contributes to this matrix element, with the result 



(88) 



The elements of the 1 x npair matrix A 



are given by 




^q+P/2^-q+P/2 



(89) 




The only contribution to this matrix element comes from through the term hi2 , giving 




^ /ii,2(q + P/2,-q 




+ P/2) 



(90) 
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Using the expression for hi2 term show in Eq. ()77|) we can write ()9Up in terms of the coherence 
factors (|H3j). (jHBj) obtaining 



A'\q; P) = , ^ {Co(+)Ci(P, + Co(-)Ci(P, +)V^(-) 

^1 + Ko 



-Co{P)C,{+,-)V{P)} (91) 

where we used the notation ± = ±q + P/2 and ±' = ±q' + P/2. The elements of the 
npairs x 1 matrix A^^ is the hermitian conjugate of the matrix A^"^. The elements of 
npairs x npairs matrix ^^^(q', q; P) are given by 

A'\c^,q;P) = ($ I , H, ] 1 .1,) , (92) 



Both /i2, through the term /ii i, and through the term /i2,2 contribute to this matrix 
element, leading to 



A'^ic^, q; P) = [^(-q + P/2) + ^(q + P/2)]<5q,q, + , ^ ^ I 

^(l+5q,0)(l+5q',0) 



/i2,2(q + p/2, -q + P/2, q' - q) + /i2,2(q' + P/2, -q' + P/2, q - q') + 

/i2,2(-q + P/2, q + P/2, q' + q) + /i2,2(-q' + P/2, q' + P/2, q + q')}. (93) 

Using the expression for /i2,2 calculated in the previous appendix we get 



A'^ic^, q; P) = [ivi-] + u;i+)]5^,^, + { 

:i+'5q,0)(l+5q',0) 



C2{+', -')C^2(+, -)V{P) + +')Ci(-, -')V{\ci - q'l) + 
C^i(+',-)Ci(+,-')V^(|q + q1)} (94) 



The matrix B 

For the matrix B we start with B^^ 



fiii(P) = ($|[r7P,i/,r/_p]|$), (95) 
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By inspection, we see that there are no terms in the Hamiltonian that contributes to this 
matrix element, hence 

B^\P) = 0. (96) 
Next we consider the matrix elements of B^^ given by 

^^^(q;P) = ($|[r;p,i/,^?=S=^a^]|$), (97) 

In this case the only contribution comes from through the term h^ Q 
i3^^(q;P) = - [/^3,o(q + P/2,-q+P/2) + 

/i3,o(-q + P/2, -P) + /i3,o(q + P/2, -P)] (98) 
which can be written in terms of the coherence factors as 

^''(q;P) = - . V {Co{+)c,{p,-)v{+) 

Co(-)C2(P, +)V{-) + Co(P)C2(+, -)V{P)]} (99) 
The matrix is the transpose of B^^. 

For S^^(q', q; P) we have 

S^^(q', q; P) = {^\[!I^±lI^Z^±l/l , H, !k±Z|2=^|$), (100) 

The only non vanishing contribution to this matrix element comes from /i4 through the term 
^4,0, giving: 

B'\c^, q; P) = -/io,4(-q + P/2, -q' - P/2, q - q') 
-ho4q + P/2, -q' - P/2, -q' - q) 
-/io,4(-q + P/2, q + P/2, -q' + q) 
-hoA<l + P/2, -q + P/2, -q - q') 
-V4(-q+P/2,q'-P/2,q' + q) 
-Mq + P/2,q'-P/2,q'-q) 
-V4(-q+P/2,q + P/2,q + q') 
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-/io,4(q + P/2,-q + P/2,q'-q) 

-hoM - P/2, -q + P/2, P) - hoM - P/2, q + P/2, P) 
- V4(-q' - P/2, -q + P/2, P) - /io,4(-q' - P/2, q + P/2, P) 

q',P ^ -P (101) 

which can be written, in terms of the coherence factors as 



S^^(q', q; P) = - , ^ ^ A C,[+, +']C,[-, -'][y(|q + q'|) 

V(l + M(l + <^q',0) 

+y(|q - q'l)] + C,[+, -]C,[+', -'MP)} (102) 
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